The QCD Critical End Point Under Strong Magnetic Fields 
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We use the three-flavor Nambu-Jona-Lasinio model, which includes strangeness and quark physi- 
cal masses in the mean field approximation, to investigate the influence of strong magnetic fields on 
the QCD phase diagram covering the whole T — \i plane. It is found that the size of the first order 
transition line increases as the field strength increases so that a larger coexistence region for hadronic 
and quark matter should be expected for strong magnetic fields. The location of the critical end 
point is also affected by the presence of magnetic fields which invariably increase the temperature 
value at which the first order line terminates. On the other hand, the critical end point chemical 
potential value displays a subtle oscillation around the B = value for magnetic fields within the 
10 17 - 10 20 G range. These findings may have non trivial consequences for the physics of magnetars 
and heavy ion collisions. 
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Magnetic catalysis is an interesting phenomenon with 
direct consequences on the QCD chiral symmetry break- 
ing (CSB) mechanism. At vanishing T and fi one ob- 
serves, within effective quark models, that the presence 
of a magnetic field (B) stabilizes the chirally asymmet- 
ric vacuum by antialigning the helicities of a quark-anti- 
quark pair (qq) which is then bound by the strong interac- 
tion. This contrasts with the behavior of the BCS ground 
state of an ordinary superconductor where B favors the 
alignment of the spins for an electron pair opposing pair 
formation. As a consequence of magnetic catalysis, the 
quark condensate, and hence the quark effective mass, as- 
sumes higher values as B increases. Predictions based on 
the Nambu-Jona-Lasinio model (NJL) [l| suggest that 
the effective quark mass raises by approximately 30 to 
40% when B increases from 10 17 to 3 x 10 19 Gauss, i.e., 
by about two orders of magnitude @, Q • ^ is important 
to stress that the effective masses at vanishing B coin- 
cide with the ones obtained for B = 10 17 G or lower. 
The fact that magnetic fields induce magnetic catalysis, 
enhancing CSB, naturally leads to the question of how 
these fields would influence chiral phase transitions at fi- 
nite temperatures and/or chemical potentials. Of course 
this is more than merely an academic question since one 
can immediately point out at least two realistic scenarios 
whose physics are influenced by the behavior of strongly 
interacting matter under intense magnetic fields. The 
first refers to the high temperature and low chemical po- 
tential regime relevant to non central heavy ion collisions 
where huge fields, of the order \eB\ > m 2 ~ 10 18 G, can 
be created by heavy ion currents due to the spectator nu- 
cleons. The field intensity depends on the centrality and 
beam momentum so that \eB\ w 5m 2 can be reached at 
RHIC while \eB\ « 15m 2 can be reached at the LHC. 
Although time dependent and short-lived 3 these fields 
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may have a non negligible effect on the transition 
Strong magnetic fields in heavy ion collisions at beam 
energies ranging from 0.2 to 2 GeV/nucleon, which cover 
most of the accelerators in the world, were considered in 
0, where the authors conclude that particles with low 
masses and their differential ratios (as n~ for exam- 
ple) can be used as possible probes of the density depen- 
dence of the symmetry energy, another topic of interest in 
the recent literature. The second scenario refers to the 
cold baryon-dense matter which forms compact stellar 
objects such as magnetars. While ordinary neutron stars 
bear a magnetic field of the order of 10 9 — 10 12 G, mag- 
netars, believed to be the source of intense gamma and 
X rays, bear fields of the order of 10 13 — 10 15 G at their 
surface reaching values up to 3 to 5 times of magnitude 
greater in their center 7]. As far as phase transitions 
are concerned, these two physically appealing situations 
are located at the high-T-low-/x (relativistic heavy ion 
collisions) and low-T-high-/x (neutron stars) extremes of 
the QCD phase diagram while the intermediate region 
lies within the capabilities of experiments such as the 
low energy scan in HIC at RHIC, FAIR at GSI, NIC A 
at JINR, and J-PARC at JAERI where, eventually, high 
magnetic fields will also be reached. 

The determination of the QCD phase diagram, even 
at vanishing B, is still a matter of great theoretical and 
experimental activities. In this case, powerful lattice sim- 
ulations have established that at vanishing baryon densi- 
ties there is no true phase transition from hadronic mat- 
ter to a quark gluon-plasma but rather a very rapid raise 
in the energy density signaling a crossover characterized 
by a pseudocritical temperature, T pc , which is expected 
to be within the 150 — 200 MeV range, with systematic er- 
rors [5j . The situation is less clear for the finite chemical 
potential region since, so far, there is no reliable infor- 
mation available from lattice QCD evaluations. Never- 
theless, most finite /j, lattice extrapolations for the \i = 
Columbia plot indicate that the critical first order surface 
(on the m Ui d — m s — /i plane) will hit the physical cur- 
rent mass values at some finite fi thereby characterizing 
a critical end point (CEP) @. 

Therefore, the use of effective models for QCD consti- 
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tutes the pragmatic approach to treat the finite chemi- 
cal potential region, especially at very low temperatures. 
Within this spirit, the three-flavor version of the NJL has 
been widely used in investigations which aim to analyze 
the phase structure of chiral transitions. Using a stan- 
dard parametrization, at B — 0, this model predicts a 
phase diagram [9] which is in accordance with what is the 
most general current belief. Namely, the model predicts 
that at vanishing chemical potential a crossover takes 
place with a T pc which is in good agreement with the 
lattice results. At zero temperature, a first order transi- 
tion occurs for /i ~ Mb /3, where Mb is the baryon mass, 
so that a first order transition line emerges at smaller fi 
as T departs from zero. Since at \x = a crossover is 
predicted, this line cannot cover the whole T — \i plane 
and will terminate at a CEP located at intermediate val- 
ues of T and //. Considering this picture as the B = 
benchmark one can further analyze how a magnetic field 
will influence the expected QCD chiral transition. So far, 
most model applications carried out with this aim were 
performed at [i = with the two-flavor linear a model 
(LSM) [HI as well as with the two-flavor NJL model [ll[ 
predicting that the crossover takes place at higher T pc 
when B ^ 0. Another interesting outcome from these 
investigations shows that the chiral and deconfining lines 
can split at finite B when the models also include the 
Polyakov loop. On the lattice side, a two flavor simula- 
tion, which used quark current masses corresponding to 
pion masses on the range 200-400 MeV has shown that 
the deconfincmcnt and chiral symmetry raise by just a 
few percent even for very high fields and have not ob- 
served any splitting However, very recently another 
lattice simulation has considered 2+1 flavors and phys- 
ical values for the quark current masses predicting that 
the pseudo critical temperature should decrease with in- 
creasing B [j| . This also constitutes an interesting result 
if one recalls that the running of the strong coupling in 
the presence of B may turn the magnetic catalysis effect 
around, as suggested in Ref. [l3j . 

Despite the progress made at vanishing baryon den- 
sity very little has been done to determine the influence 
of magnetic fields on the whole T — u plane. One such 
investigation was carried out in Ref. [14| using the stan- 
dard two flavor LSM but when the vacuum is included, 
as it should be [l(| , the whole T — \x plane is dominated 
by the crossover just as it happens at B = 0. On the 
contrary, as shown in (l5| . the standard two flavor ver- 
sion of the NJL model in the chiral limit does not share 
this property allowing for the observation of how a mag- 
netic field influences the location of the tricritical point. 
In Ref. [l5j a Fock-Schwinger proper time method was 
used in the chiral limit. The authors have observed a cu- 
rious pattern between magnetic catalysis and the Haas- 
van Alphen effect, which leads to oscillations in the first 
order transition line. The complexity of the phase dia- 
gram obtained with several tricritical points is probably 
due to the zero mass quark limit and further investiga- 
tion with realistic current quark masses is needed. As 



emphasized in (l6| . strange quarks play a crucial role in 
shaping the phase diagram of QCD since their mass, m s , 
controls the nature of the chiral and deconfinement tran- 
sition. At the same time, this mass has also an important 
effect on the stability limit of neutron stars and on the 
possible existence of a quark core in collapsed stars. For 
these objects, the establishment of the critical first or- 
der line is of utmost importance once the possibility for 
hybrid and quark star formation depends on the quark 
matter nucleation process. The determination of the first 
order line, at B ^ 0, may also be relevant in the frame- 
work of studies related to thermal nucleation of quark 
droplets in hadronic matter found in protoneutron stars 
whose temperatures are of the order 10-20 MeV. 

In this Letter we report on the main results obtained by 
scanning the whole T — \i plane paying special attention 
to the CEP and to the first order coexistence region. 
We do not consider the Polyakov loop since its effects 
are more pronounced at the high temperature and low 
chemical potential region was shown 

that, although it changes considerably the values of the 
temperature and chemical potential related to the CEP, 
the qualitative behavior of the phase transition was not 
modified. 

By mapping the first order line on the T — fi plane 
we can automatically predict the location of the CEP for 
different values of B. At vanishing \i we reproduce the re- 
sults known from the other model applications (like those 
in Refs. fTol. flTjj ^ predicting that the B = crossover will 
take place at a higher T pc for B ^ 0. The new results re- 
ported here emerge at the finite baryonic density domain 
where we observe that the presence of a magnetic field 
favors the presence of latent heat at higher temperatures, 
invariably increasing the size of the first order transition 
region so that the CEP temperature, Tcep(-B), increases 
with B. We also show that the critical chemical poten- 
tial, /x cr it(-B), which determines the low temperature first 
order transition displays a non trivial behavior moving 
around the B = value for magnetic fields in the range 
10 17 - 10 20 G. To the best of our knowledge this is the 
first time such predictions are made within a more re- 
alistic approach, which takes strangeness and physical 
masses into account. 

In order to consider (three flavor) quark matter sub- 
ject to strong magnetic fields we introduce the following 
Lagrangian density 

C = Cf-\F liV F^, (1) 

where the quark sector is described by the su(3) version 
of the Nambu-Jona-Lasinio model which includes scalar- 
pseudoscalar and the t'Hooft six fermion interaction that 
models the axial U{1)a symmetry breaking [l8j : 

Cf=^ f [j fl (id" - qfA*) - rh c ]ipf + C sym + C det , (2) 
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FIG. 1: a) Critical chemical potential in terms of the magnetic field for different values of the temperature; b) Phase diagram 
on the T — fj, plane. The thick curves represent first order transition lines which terminate at a critical end point identified 
with a full dot and the thin lines represent a crossover. 



with C S ym and C de t given by [ijj]: 
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C sym = GY / [^f^f) 2 + (^fil 5 Xa^f) 2 ] , (3) 
a=0 

C det = ~K {det/ + 75)^/] + det/ [1/7 (1 - 75)^/] } 

where ipf = (u, d, s) T represents a quark field with three 
flavors, rh c = diagy (m u , rrid, m s ) is the corresponding 
(current) mass matrix while qf represents the quark elec- 
tric charge, Ao = \/2/3I where I is the unit matrice 
in the three flavor space, and < A a < 8 denote the 
Gell-Mann matrices. The model is not renormalizable, 
and as a regularization scheme for the divergent ultra- 
violet integrals we use a sharp cut-off, A, in 3-momentum 
space. The parameters of the model, A, the coupling 
constants G and K and the current quark masses to 9 
and m 9 are determined by fitting /„., , ttik and 
rrifji to their empirical values. Throughout this paper 
we consider A = 631.4 MeV , m u = rrid = 5.5 MeV, 
m s = 135.7 MeV, GA 2 = 1.835 and KA 5 = 9.29 as in 

As usual, and F M „ = d^A v — d u A^ are used to 
account for the external magnetic field. We consider 
a static and constant magnetic field in the z direction, 
A^ = S^xiB. 

The thermodynamical potential for the three flavor 
quark sector, Qf is written as 

n f {T,fi,B) = 2G tf-^KfafacPk 

i=u, d, s 

+ (n» f ac + «™ a9 + nj ed ^j , (5) 

where the vacuum (f2™ c ), the magnetic (0™ aff ), the 
medium contributions (£l™ ed ) and the condensates <f>i 
have been evaluated with great detail in Ref. [2, 3]. The 
effective quark masses are solutions of the gap equations 



obtained from the minimization of the thermodynamic 
potential with respect to the effective masses, 

M i = m i -AG<t> i + 2K<t> j <t>k> (6) 

with (i,j,k) being any permutation of (u,d,s). In this 
work we consider fi = (A u = Mrf = Ms • 

The first order phase transition line is obtained by en- 
forcing the Gibbs conditions while the cross over region 
is determined by the maximum of —d(j> m /dT [m = u, d). 

In Fig. QJi we plot the critical chemical potential as a 
function of the magnetic field for different low tempera- 
ture values where the first order transitions take place. 
For a fixed magnetic field, the critical chemical potential 
(Merit) is approximately constant up to magnetic fields of 
the order of 10 19 G, where it reduces drastically reaching 
a minimum value around 3 x 10 19 G and then it increases 
again. This occurs precisely when just one Landau level 
is filled. One can clearly sec the effects of the filling of 
the Landau levels at low temperatures, where the van 
Alphcn oscillations are seen. 

The effects of the behavior observed in Fig. [IJl are re- 
flected in the curves shown in Fig. [TJj, where a T — fi 
phase diagram is shown for different values of B. As the 
magnetic field increases, the line of the first order phase 
transition also increases. The line moves towards a lower 
chemical potential and then back to higher chemical po- 
tentials for increasing values of B. The highest point in 
each first order transition curve is the CEP, indicated 
with a full dot. The curve obtained for B = 10 17 G 
coincides exactly with the results presented in [9( for a 
zero magnetic field. The pattern observed in Fig. [TJd 
can be seen with more detail in Fig. [2] where the depen- 
dence of hcep and Tcep, which characterize the CEP, 
are shown as a function of the magnetic field B. Fig. 
[5^ shows that hcep essentially decreases with B, with a 
lower decrease for fields stronger than 10 19 G. Just be- 
low 10 20 G, hcep starts to increase. Fig. [5b represents 
fiCEP versus Tcep, each point corresponding to a differ- 
ent value of the magnetic field: the smallest value of the 



4 




60 80 100 120 140 160 



T CE p(MeV) 



FIG. 2: Critical end point chemical potential \icep in terms of a) the magnetic field, and b) the critical end point temperature 
Tcep; c) Tcep as a function of B. 



critical chemical potential, 258.5 MeV, corresponding to 
a T CE p ^ 137 MeV occurs for B = 8 x 10 19 . On the 
other hand, Tcep always increases with the increase of 
B, mainly above 10 19 G. 

In summary, our results show that, at fi — 0, T pc in- 
creases with B in agreement with most model calcula- 
tions [13, [H| an d an early lattice simulation [l2[ which is 
not surprising since we are dealing with the NJL within 
the MFA. On the other hand, the main contribution of 
our work concerns the critical end point and the first or- 
der coexistence region for which magnetic effects become 
noticeable for fields greater than 10 17 G. In this case, 
one observes an increase on the size of the first order 
transition line with Tcep moving to higher values as B 
increases. This result can be of interest for the physics 
of heavy ion collisions at intermediate energies. We have 
also observed that, due to the van Alphen effect, the first 
order transition lines display a non trivial behavior os- 
cillating around the B = line. This implies that, as 



B increases, criticality will be achieved at chemical po- 
tential values which can be greater or smaller than the 
B = value, /i cr it ~ 360 MeV. These findings may have 
direct consequences regarding the dynamics of the first 
order phase transition which influences the quark matter 
nucleation process. Although we have attempted to per- 
form a realistic application, by considering three quark 
flavor and physical mass values, it is important to recall 
that at temperatures higher than 100 MeV the effects of 
the Polyakov loop may be important 21]. The Polyakov 
loop field is going to be more carefully considered in a 
forthcoming paper but we belive that these seminal qual- 
itative results are not going to be altered. 
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